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To test a possible relation between the topological entropy and the Arnold complexity, and to 
provide a non trivial example of a rational dynamical zeta function, we introduce a two-parameter 
family of two-dimensional discrete rational mappings. The generating functions of the number of 
fixed points, and of the degree of the successive iterates, are both considered. We conjecture rational 
expressions with integer coefficients for these two generating functions, and a rational expression 
for the dynamical zeta function. We then deduce algebraic values for the complexity growth and 
for the exponential of the topological entropy. Moreover, these two numbers happen to be equal for 
all the values of the parameters. These conjectures are supported by a semi-numerical method we 
explain. This method also localizes the integrable cases. 
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Chaos is associated to extreme complexity and "unpre- 
dictability" . Very few exact results are therefore known, 
however several exponents have been introduced to pro- 
vide measures of the complexity and further, to classify 
chaotic systems |^J^. The most popular are the Lya- 
pounov exponents, which have a clear and intuitive in- 
terpretation, but require the choice of a metric and of 
an invariant measure. They can vary considerably un- 
der a very tiny change of the parameters of the dynam- 
ical system (see |^ p237): probing fine details, they are 
not universal. The same remarks also apply to the met- 
ric entropy By contrast there also exist exponents 
which do not involve any assumption on the phase space. 
These exponents, obviously, give a less detailed descrip- 
tion of the system, but are more universal. They provide 
a mean for a classification of dynamical systems. The 
topological entropy ||] log h and A the Arnold complex- 
ity Q are two examples. The topological entropy probes 
the growth of the number of stable cycles as a function 
of their length, and the Arnold complexity probes the 
growth of the number of intersections of a line with its 
successive iterate. These two notions give general infor- 
mations not sensitive to specific details. From an intu- 
itive point of view one can understand the importance of 
the topological entropy, since the asymptotic behavior of 
a dynamical system heavily depends on its fixed points 



and stable cycles: an initial point is often very close to 
many basins of attraction, which results in a chaotic mo- 
tion. Precisely we introduce the fixed points generating 
function H(t) = hn • where ft,„ is the number of 
real, or complex, fixed points of the n-th power A:" of the 
mapping. The same information is also coded in the so- 
called dynamical zeta function C,{t) introduced by 
Artin and Mazur |^ and related to the generating func- 
tion H{t) by H{t) = t^ln(C(i))- Both functions only 
depend on the number of fixed points, and not on their 
particular properties or localization: functions H{t) and 
t^(t) are invariant under topological conjugacy (see Smale 
JlOt for this notion). They do not depend on a specific 
choice of variable, h, the exponential of the topological 
entropy characterizes how the coefficients /i„ grow with 
n: hn ~ h", so that h is the modulus of the inverse of 
the smallest modulus pole of H(t), if rational. For a lin- 
ear dynamic on a torus, the cat map, the exponential of 
the topological entropy has been calculated and found 
algebraic, i.e. solution of polynomial with integer coef- 
ficients. We are not aware of other non trivial dynami- 
cal systems where an algebraic value for the exponential 
of the topological entropy has been calculated. In this 
letter we will provide such an example of a discrete dy- 
namical system with a rational dynamical zeta function, 
and, consequently, algebraic value for the exponential of 
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the topological entropy. This result should not be con- 
sidered like a mathematical curiosity : it is similar to 
the rationality of critical exponents in conformal theory. 
This algebraicity is the signing of deeper "rigid" struc- 
tures (like Feigenbaum cascades |jll| are). 

In the context of rational mappings it is easy to see 
that the Arnold complexity can be replaced by the com- 
plexity growth of the successive iterations. To define the 
complexity growth A we introduce the complexity gener- 
ating function G{t) — J^n ^" ' ^" where d„ is the degree 
of any of the numerator, or denominator, of the compo- 
nents of the successive iterates of the rational mapping 
under consideration. When common polynomials factor- 
ize in the numerators and denominators, the coefficients 
dn grow slower than expected. We stress that this defi- 
nition only apply to rational mappings. The complexity 
growth A characterizes how coefficients dn grow with n: 
dn ^ A". Like h, the exponential of the topological en- 
tropy, complexity A is the modulus of the inverse of the 
smallest modulus pole of G{t). In this letter we claim 
that the complexity growth A and h, the exponential of 
the topological entropy, are equal 



h = X 
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This will be tested successfully for a particular class of 
mappings, for which both generating functions are con- 
jectured to be rational and, consequently, the complexity 
growth and the exponential of the topological entropy are 
algebraic. We will also give an effective semi-numerical 
method to compute these two characteristic numbers. 

Let us introduce the discrete rational mapping k^^e 
which associates {un+i,Vn+i) to {un,Vn) 

Un+l = 1- Un+ Un/Vn (2) 
Vn+1 = e + Vn ~ Vn/Un + a ■ (1 - Un + M„/w„) 

This mapping originates from the study of the symme- 
tries of models of lattice statistical mechanics De- 
pending on the actual values of the parameters a and 
e, the mapping can have completely different behaviors. 
For example, for e = and whatever a, as well as for 
a — Q and e — —1, , 1/2, 1/3 or 1, the mapping is 
integrable, whereas for all other values it is not ||l^. A 
simple calculation shows that ka^^ is invertible and that 
its inverse is also rational. This property of birationality 
is of importance in our study. We have formally calcu- 
lated the successive powers of ka^c for arbitrary a and e, 
from which we propose 



Ga,e(t) 



1 - t - 2i2 _ i3 



(3) 



The expansion of the conjectured expression Eq. co- 
incides with our results up to the largest power n = 7 we 
were able to compute. Another rational expression with 
the same denominator is also obtained if one uses a ma- 
tricial representation of the mapping Eq. [|l6|. The 



expression Eq. (||) of Ga,t yields A ~ 2.14789. To sup- 
port this conjecture, we have devised a semi-numerical 
method to estimate complexity A. It consists in iterating 
fca^e over the field of rationals. During the first steps, 
some 'accidental' cancellations between numerators and 
denominators can arise, but after this transient regime, 
the numerators and denominators get extremely large, 
and cancellations are only due to formal simplifications. 
We then determine how the magnitude of the four numer- 
ators, or denominators, grows with n. With this method 
it is possible to raise ka.e to the 15-th power, and more- 
over it is easy to scan a large number of values of the pa- 
rameters a and e. The calculations are performed with an 
infinite precision C-library ||l7|] . Obviously, this method 
works only for rational mappings. On Fig. [l], one clearly 
sees that, for most of the values of e, the complexity A is 
extremely close to the expected value. We call "specific" 
the values of the parameters for which the complexity 
A is different from 2.14789, they will be discussed later. 
We also have formally computed for arbitrary a and e 
the fixed points of the powers of k^^t using, once again, 
the rationality of the mapping. This gives 



iJa,e(i) = 2t + 2t^ 
+A7t^ + 
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(4) 
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From this expression we propose the rational expression 
for the generating function of the number of fixed points 



(2 + 3^2-^ t^)-t 

(l-t2).(l_i_2t2_i3) 



or, equivalently, the dynamical zeta function reads 

{i + m-e) 



{l-t~2t^-fi) 



(5) 
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Note that the total number of fixed points of fc^ ^ does not 
depend on the actual generic values of a and e, however 
the number of real fixed points is extremely dependent 
on these two parameters. Let us also mention a local 
area preserving property: the determinant of the Jaco- 
bian of the n-th power of ka.e evaluated at each fixed 
points of k^ e is equal to one. The "visual" complexity of 
the phase diagram takes its origin in the real fixed points, 
and therefore varies considerably with a and e [p^ . One 
sees that the two polynomials giving exponents A and 
h are the same, and consequently we have the equality 
h — X. This equality holds for generic values of the pa- 
rameters, however, as shown on Fig. ^, there exist specific 
values. These specific values include e — 1/3,1/2,3/5. 
It is then natural to investigate if the equality of the 
complexity growth and the exponential of the topologi- 
cal entropy also holds for the specific values. We have 
performed calculations for these values and found that 
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equality is always true. The polynomials giving the 
value of A and h are presented in table Tab. ||. Probably, 
other specific values exist, but they lead to simplifications 
occurring at very high orders, and the corresponding A 
is too close to the generic value to be distinguished from 
it with our method. 

Besides the specific values mentioned above, extra sim- 
plifications also happen for a = 0, and the complexity is 
further reduced. We hence study this special case a = 0. 
In that case, a change of variables turns /cq.e into a 
simpler mapping, 



Vn+l 
Zn+l 



Vn ■ 



+ f 



(7) 



From now on, the degrees, and the fixed points, arc those 
of fcg. Since the complexity is lower, the semi- numerical 
method presented above is more efficient and it is pos- 
sible to perform calculations beyond the 20-th power. 
The results are displayed on Fig. ^, where the existence 
of integrable values, and non generic values, is clearly 
seen. It is simple to see that, if there is no simplifi- 
cation, dn+i = dn + dn~i where d„ was introduced in 
the definition of generating function G{t). In that case 
G,{t)-2t-l = t- {G,(t) - l)+t^ -Cit). Up to the 20- 
th power there is no simplification and consequently we 
conjecture that, except for the specific values, the gen- 
erating function of the complexity growth for fcf is the 
following rational expression 



OS) = 



l + t 



1-t-P 



(8) 



The corresponding complexity growth is A ~ 1.61803, 
in excellent agreement with Fig. ||. We have studied 
the possible equality between h and A for the exam- 
ple e — 13/25 = 0.52. We have chosen this value, 
for which we present a detailed analysis, because it is 
generic. We give in table Tab. | the number of fixed 
points, as well as their properties. The corresponding 
phase portrait is very complicated and dominated by the 
real fixed points |lj] which are all saddle or elliptic. We 
note that the same properties also holds for the complex 
fixed points. The expansion of can be deduced, up 
to order eleven, from Tab. ^. This expansion is compati- 
ble with the very simple rational form for the generating 
function of the number of fixed points for kg 



Heit) 



or, equivalently, the dynamical zeta function is 

1-e 



Ut) = 



1-t-t^ 



(9) 



(10) 



As expected, the two polynomials determining the ex- 
ponential of the topological entropy and the complexity 



growth are equal, and so are A and h. Both are algebraic 
numbers. 

Coming back to Fig. |^, we now analyze, for alpha — 0, 
the specific values of epsilon. Let us recall that e = 
— 1,0, 1/3, 1/2, 1 lead to integrable mappings [|l^. This 
corresponds to a polynomial growth of complexity and of 
the number of fixed points, that is, X ^ h = I. This is 
seen on Fig. H, except for e — —1, which is out of scale 
but for which this is also true. The other specific values 
are non integrable and can be partitioned in two sets: 
{1/m; TO > 3} and {(m — l)/(m -I- 1); m> 3}. In all 
cases the polynomials giving the complexity growth and 
the exponential of the topological entropy are the same. 
These polynomials are listed in Tab. || 

In conclusion, we have given an example of two- 
parameter family of two-dimensional discrete dynamical 
system with rational dynamical zeta function and ra- 
tional degree generating function G{t). On this exam- 
ple the exponential of the topological entropy and the 
Arnold complexity have the same algebraic value. A 
semi-numerical method, applying to rational transfor- 
mations only, has been given, which enables to calcu- 
late the complexity growth, and to localize possible inte- 
grable points. In fact, and this will be detailed in forth- 
coming publications, this mapping belongs to a "huge" 
family of transformations, for which similar results are 
also obtained. This family of transformations is so large 
that (if one believes in "some" universality of dynami- 
cal systems) most of the dynamical systems should be 
very closely "approximated" by transformations having 
algebraic complexity values. 
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FIG. 1. The complexity growth A as a function of e for 
a — 27/20. e is taken of the form j/720. The arrow indicates 
the conjectured generic value. 
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FIG. 2. The complexity growth A as a function of e for 
a = 0. 

e is taken of the form j/720, the values 
a = 1/7, 1/11, 1/13, 5/7 

have also been added. The arrow indicates the conjectured 
generic value. 
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# elliptic 


1,0 


0,0 


0,0 


0,0 


1,0 


0,0 


1,0 


0,0 


1,0 


1,3 


1,0 


# saddle real 


0,0 


0,0 


1,0 


1,0 


1,0 


0,2 


1,2 


1,4 


3,4 


1,6 


5,12 


# total fixed points 


1 


1 


4 


5 


11 


16 


29 


44 


76 


121 


199 



TABLE 1. Number of real (first number) and complex (second number) fixed points of ^13/25- n-cyc\e means cycle of 
minimum length n. 
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a generic 


1 -t-t'' -2f' - t"" - t" 


1 -t -t" -2f' -t"" - 2r' -t" -t' 


generic see ([ 




(*) 


Q = 


t is n-th root of unity 


t is n-th root of unity 


i~t-e + 


1-t-i^ -t^'^+' 



TABLE IL The polynomials giving A and h in various specific cases. The symbol(*) means that a 7^ and e = (m— l)/(m+ 1) 
are not generic, however the exponents are extremely close to the generic value, preventing us to compute them reliably. The 
case a 7^ and e = 1/m is generic for m > 3. 
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